Testing Gravity with the Stacked Phase Space around Galaxy Clusters 
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In General Relativity, the average velocity field of dark matter around galaxy clusters is uniquely 
determined by the mass profile. The latter can be measured through weak lensing. We propose a new 
method of measuring the velocity field (phase space density) by stacking redshifts of surrounding 
galaxies from a spectroscopic sample. In combination with lensing, this yields a direct test of 
gravity on scales of 1-30 Mpc. Using N-body simulations, we show that this method can improve 
upon current constraints by several orders of magnitude when applied to upcoming imaging and 
redshift surveys. 
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The accelerated expansion of the Universe is the most 
tantalizing problem in modern cosmology. Within Ein- 
stein's General Relativity (GR), the cosmic acceleration 
can be explained by introducing a mysterious smooth 
component, dark energy. However, it can also be inter- 
preted as signature of the breakdown of GR on cosmo- 
logical scales. Many on-going and upcoming wide-area 
galaxy surveys aim at testing dark energy and modified 
gravity scenarios as the origin of cosmic acceleration. 

Cosmological probes of gravity are based on recon- 
structing the perturbations in the space-time metric and 
their relation to matter 0, Q . Weak gravitational lens- 
ing provides a clean measurement of the lensing poten- 
tial, while the timelike potential can be probed through 
the modulations in redshift caused by peculiar velocities 
of galaxies. In this Letter, we propose a new method 
of testing gravity at intermediate scales (1 — 30 Mpc) by 
measuring a projection of the phasespace around massive 
galaxy clusters. If GR is valid, the phasespace around 
the sampled clusters is uniquely determined by the mass 
density profile. In other words, comparing the measured 
mass density and velocity profiles allows for a model- 
independent test of Einstein gravity. The scales probed 
are complementary to and potentially provide more in- 
formation than the linear regime studied in most pre- 
vious studies 0, 0, or the small scales considered in 
[1, 0, 0| ■ Moreover, this test is a generic probe of grav- 
ity: adding other, non-standard ingredients such as mas- 
sive neutrinos or primordial non-Gaussianity, for exam- 
ple, will likely have a negligible impact on the relation 
between mass density and velocity profiles. This is not 
the case for other commonly considered probes of gravity, 
such as the matter power spectrum or cluster abundance. 
The main challenge lies in modeling the observables on 
these scales. We will demonstrate the feasibility of our 
method by using N-body simulations for Einstein and 
modified gravity models. 

Methodology: Consider a sample of galaxy clusters 
(with accurate redshifts) in a cosmological volume cov- 
ered by a spectroscopic galaxy survey. We can then con- 
struct the two-dimensional distribution of galaxy-cluster 
pairs in terms of the transverse distance r p , and the rel- 



ative line-of-sight velocity «i os . More precisely, we have 
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where A9 gc is the angular separation of galaxy and clus- 
ter, dji is the comoving angular diameter distance, c is the 
speed of light, and z g , z c denote the galaxy and cluster 
redshifts, respectively. The average phase space distribu- 
tion is estimated by stacking all cluster-galaxy pairs. 

The lower panel of Fig. [1] shows this distribution mea- 
sured in the Einstein-gravity N-body simulations of [Tl| 
around halos with masses M > 10 14 M Q /h identified at 
z = 0.35, where we assumed a concordance A and cold 
dark matter cosmological model (ACDM). We use the 
output at z = 0.35 of 20 simulations of 1.5 (Gpc/h) 3 vol- 
ume each. We defined halos using the friends-of-friends 
finder algorithm with linking length 0.2 times the mean 
particle separation, and assigned center-of-mass positions 
and velocities using the member particles. To mimic a 
galaxy redshift survey, we select secondary halos in the 
mass range 3 x 10 13 < M s < 10 14 M Q //i as galaxies in 
a cube of side length 40 Mpc/h centered on each pri- 
mary halo. In real galaxy surveys, such a selection in 
real space is not possible of course; we will return to this 
point below. By stacking over many clusters and binning 
in r p , we average over triaxial or irregular density pro- 
files, yielding a distribution which is only a function of 
r p and |v los |. 

The lower panel of Fig. [1] clearly shows two distinct 
regimes; at small radii r p < 2Mpc//i, iso-density con- 
tours are closed, while on larger scales the contours be- 
come open, reflecting the ongoing infall onto the mas- 
sive halos. The boundary between these two regimes has 
been used in the caustic method [12]. As shown in a 
forthcoming paper, we can construct an accurate model 
of the Vios — T p distribution of dark matter halos through 
a combination of N-body simulations and analytical the- 
ory. The upper panel shows the analytical prediction 
for the RMS dispersion of v\ os as a function of r p (esti- 
mated through the standard sample RMS). The model 
prediction is in good agreement with the simulation re- 
sult, within the statistical errors of the simulation mea- 
surements, which were measured from 20 simulation re- 
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FIG. 1: Lower panel: The v\ os — r p phase space distribu- 
tion (in logarithmic scale) as measured using halo catalogs 
constructed from N-body simulations in ACDM; we consid- 
ered primary halos ( "clusters" ) with masses > W 14 M Q /h and 
secondary halos ("galaxies") in the range 3 x 10 13 < M < 
10 14 M@/h. Upper panel: The dispersion of the line-of-sight 
velocity distribution cr„ l0B as function of r p . The data points 
with error bars are computed from the simulation results in 
the lower panel, while the solid curve is our analytical model 
prediction. The error bars are scaled to mimic the measure- 
ment accuracies for a spectroscopic survey of 2000 sq. degrees 
over 0.2 < z < 0.4. 



alizations so as to mimic the measurement accuracies for 
a survey of 2000 sq. degrees coverage and with rcdshift 
range 0.2 < z < 0.4 (see below for more details). How- 
ever, there are two complications that in reality need to 
be taken into account: the contribution to vi os from the 
cosmological redshift, which is given by HAr\ os , where 
An os is the line-of-sight separation between the galaxy 
and the cluster; and the contribution from motions of 
galaxies within their parent halos. The Hubble flow con- 
tribution can be modeled if the real-space cluster-galaxy 
correlation function on scales of interest is known. In 
practice, we can only measure the redshift- space correla- 
tion function, which in turn receives contributions from 
the velocities. One approach is to construct a joint model 
of the vi os — r p phase space and the redshift-space correla- 
tion function. Another possibility is to measure stacked 
weak lensing around the galaxies, which yields the real- 
space galaxy-matter correlation function. Combining the 
cluster-matter and galaxy-matter correlation can be used 
to infer the galaxy-cluster correlation function. 

Another effect which has to be included is the motion 
of galaxies relative to the center-of-mass of their parent 
halos. In order to include this contribution, we need to 
know the distribution of relative velocities as well as ra- 
dial offsets relative to their halos. If the galaxies are 
dynamically relaxed within the halos, these distributions 
are related by the virial theorem. Stacked weak lensing 
measured for the galaxy sample yields the mean parent 
halo mass as well as giving clues to the distribution of 
radial offsets. This can be used to constrain the galaxy 
motions within halos [l3j . 

An accurate measurement of the phase space distribu- 
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FIG. 2: Upper panel: Ratio of the velocity dispersion o~ v 
along the line of sight measured around halos with M300 > 
10 14 Mq//i in f(R) simulations to that measured around ha- 
los of the same mass in ACDM simulations. Lower panel: 
Ratio of the projected density profiles of the same halos in 
f(R) and ACDM simulations. This is approximately what 
stacked lensing would measure. The error bars in both pan- 
els are estimated from the simulations, as in Fig. [T] for the 
spectroscopic and imaging surveys of 2000 sq. degrees. 
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FIG. 3: Same as Fig. H but for DGP models (see text). 



tion requires a redshift survey of galaxies which covers a 
sufficiently large cosmological volume to achieve a given 
number of cluster-galaxy pairs. One further advantage of 
the stacking procedure and of considering scales of sev- 
eral Mpc is that we do not necessarily require a high 
number density of spectroscopic galaxies. In contrast, 
deep dedicated observations would be needed if one were 
to determine the velocity dispersion of individual clus- 
ters. The stacked weak lensing measurement requires an 
adequately deep imaging survey so that images of back- 
ground galaxies are well resolved. 

Results: We now turn to the signatures of modified 
gravity in the t>i os — r p phase space. We begin with 
the modified action f(R) model (see (Tij l and references 
therein), specifically the one of [15( with n = 1. The 
model can be parametrized by the amplitude fuo of the 
scalar degree of freedom fu = df/dR today, with /ro = 
being equivalent to ACDM. For the values considered 
here (\f m \ = 10~ 4 - 10~ 6 ), the expansion history is in- 
distinguishable from ACDM. Current best cosmological 
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constraints arc a \f R0 \ <a fewxlO -4 [H, G3- In f{R) 
gravity, gravitational forces are enhanced by a factor of 
4/3 within the redshift-dcpcndcnt Compton wavelength 
of the field. In addition, this model incorporates the 
chameleon mechanism which restores GR in high-density 
environments [lH EH- The upper panel of Fig. [2] shows 
the dispersion a v of the line-of-sight velocity distribution 
in bins of r p measured around halos above 10 14 M /h in 
f[R) N-body simulations [HI, [2(|, relative to that mea- 
sured in ACDM simulations around halos above the same 
mass threshold [29]. Due to the limited volume and reso- 
lution of the modified gravity simulations, we performed 
this measurement for dark matter particles. Note that 
the enhancements in a v can become significantly larger 
than the effect on the virial velocities, which are enhanced 
by up to a factor of « 1.15 in f(R) 0]. 

In case of chameleon theories such as f(R), if the spec- 
troscopic galaxies are screened, we expect the enhance- 
ment of velocities to be suppressed . Secondary halos 
with AI 3Q0 > 3 x 10 13 M Q /h identified in the f[R) sim- 
ulations indicate a somewhat suppressed effect on the 
velocity dispersion for f R Q < 10 -5 , although the error 
bars are large. This is consistent with the mass thresh- 
olds ~ 10 14 M //i and below for the chameleon mecha- 
nism for these field values jij. On the other hand, the 
chameleon-screening of the clusters only affects the phase 
space at separations of order the virial radius of the clus- 
ter halos, i.e. a few Mpc or less. This can be seen for 
the cases of f RQ = 10~ 5 , 10 -6 in Fig. [2j for which the 
primary halos are screened in the simulations. 

Fig. El shows the same measurement in N-body sim- 
ulations of the Dvali-Gabadadze-Porrati (DGP) type 
braneworld models pll . I22I ] . We consider simulations for 
a self-accelerating DGP model without any A or dark 
energy (sDGP, [231]). and normal-branch models includ- 
ing a dark energy component (nDGP, [24j|). The dark 
energy equation of state is adjusted to yield a ACDM 
expansion history, making these nDGP models indistin- 
guishable by geometric probes 



24|. In case of sDGP, 



we compare to a GR model with an effective dark en- 
ergy yielding the same expansion history, in order to iso- 
late the modified structure growth effects. DGP mod- 
els are characterized by the cross-over scale r c , above 
which gravity transitions from 4D to 5D. On scales be- 
low r c , gravity is described by a 4D scalar-tensor theory, 
where the strength of the modified force scales with Hr c . 



In sDGP. 



1.35-Efn 



4038Mpc//i = 6118 Mpc, 



while in nDGP-1 (-2) it is taken to be 500 (3000) Mpc. 
As expected, we see that sDGP yields smaller veloc- 
ities than GR, since gravity is weakened in the self- 
accelerating branch. Conversely, normal-branch models 
yield higher velocities. We find no indication of a sup- 
pression of the effect when considering secondary halos 
[M s > 3 x 10 13 Mq/K) instead of dark matter, consistent 
with the fact that the Vainshtcin screening mechanism 
inherent in these braneworld scenarios does not directly 
lead to a velocity bias ||. 

Will upcoming surveys be able to detect such modi- 



fied gravity signatures in the phase space distribution? 
The statistical uncertainties in a v arise from an imper- 
fect sampling due to a finite number of the cluster-galaxy 
pairs and from cosmic variance due to a finite volume 
coverage. To make realistic forecasts, we adopt survey 
parameters that resemble the planned imagin g an d spec- 
troscopic surveys with the Subaru Telescope [l3j; we as- 
sume a survey area of 2000 square degrees, and consider 
as cluster sample halos with mass greater than 10 14 M Q //i 
and in the redshift range 0.2 < z < 0.4. The comov- 
ing volume corresponds to 0.23 (Gpc//i) 3 . We chose the 
mass range so that the massive halos allow an accurate 
measurement of the average mass profile with weak lens- 
nig [26|. We choose a cluster sample at relatively low 
rcdshifts in order to allow for a denser sampling of red- 
shifts of the secondary halos (galaxies). For the latter, 
we assume that the galaxies reside in halos with masses 
3 x 10 13 < M s < 10 14 M Q /h as in Fig.Q] and assume one 
galaxy per halo residing at the halo's center of mass. The 
mean number densities of the primary and secondary ha- 
los, found from the ACDM simulations, are 1.7 x 10 -5 
and 8.3 x 10 -5 [Mpc//i] -3 , respectively, the latter be- 
ing lower than the density of spectroscopic galaxies for 
the SDSS BOSS survey [27| or the target density for the 
Subaru survey. 

To estimate the measurement accuracies, we divide 
each of our 20 realizations of N-body simulations for 
ACDM into 27 subvolumes of 0.056 (Gpc//i) 3 at the 
output redshift z = 0.35, in order to increase the sam- 
ple size. We compute the stacked phase space distribu- 
tion [v\ os ,Tpj and RMS <r v {r p ) using pairs in each sub- 
volume. We then compute the mean and covariance of 
the ovprofiles for all radial bins over the 540 samples. 
Finally, we rescale the covariance by (Kdm/Kmrvey) 1 ■ 
The la uncertainty in each radial bin is shown in the 
upper panel of Fig. [5] and [3J The constraining power 
of the assumed galaxy survey is clearly very significant, 
over a wide range of separations. Note that the er- 
ror bars at different radial bins are highly correlated. 
To be more quantitative, we can estimate the value of 
Ax 2 = -2Aln£ between the ACDM and f(R) models, 
using the full covariance of <J v (r p ) as measured in the 
ACDM simulations. This yields Ax 2 = 218, 70 and 2.2 
for |/fio| = 10~ 4 , 10 -5 and 10~ 6 , respectively, assuming 
that the shape of the velocity profile is perfectly known. 
This shows that there is enough signal-to-noise to probe 
f[R) gravity down to field values at which the secondary 
halos become chameleon-screened and this measurement 
loses its power. Adding a log-normal scatter in mass 
of cr(ln(M)) = 0.2 in both primary and secondary ha- 
los changes the velocity profile by less than 5%. Hence 
the velocity profile appears to be robust with respect to 
uncertainties in mass estimates of the halos. 

For comparison, the lower panels of Fig. [2] and [3] show 
the ratio of the stacked spherically averaged profile of the 
mass density around the primary halos in modified grav- 
ity to that in ACDM. This quantity can be reconstructed 
from weak lensing measurements, and has been used to 
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constrain f(R) gravity [28|. Clearly, the departures in the 
mass profile are much smaller than those in the veloci- 
ties. The error bar at each radial bin shows the expected 
la measurement uncertainties for a Subaru-type imag- 
ing survey covering the same region of the sky, i.e. 2000 
square degrees. The lensing errors are determined by the 
survey area and the shot noise [2(| ; we assumed a back- 
ground galaxy density at z s > 0.6 of n g = 22 arcmin -2 
and a RMS intrinsic ellipticity of er e = 0.22. Given the 
size of the error bars relative to the modified gravity ef- 
fects, it is clear the lensing signal itself is a much less 
powerful probe of gravity than velocities. 

Discussion: In this Letter, we have investigated a 
method of using the phase space distribution around 
massive clusters to constrain modified gravity models. 
Using collisionless numerical simulations for ACDM, 
f(R) and DGP models, we demonstrated that the 
velocity dispersion as a function of transverse separation 
shows up to order unity deviations when compared to 
the profile in GR. On the other hand, the effect on the 
mass density profile, which is measurable through weak 
lensing, is much less affected by modifications to gravity. 
While we have concentrated on the second moment of the 
velocity distribution here, in principle even more infor- 
mation is contained in the higher moments. As working 
examples, we showed that a spectroscopic survey cover- 
ing an area of 2000 square degrees can in principle yield 
greatly improved constraints on f(R) and DGP models 
(see Figs. [2] and [3]). The scales probed by this method 



arc in the (weakly) nonlinear regime, and bridge the gap 
between the scales probed by redshift-space distortions 
in galaxy two-point correlations on large scales and virial 
velocities within halos on small scales. By combining 
these different methods, we can probe gravity properties 
over a wide range of scales, and have a better chance 
of capturing the signatures of the screening mecha- 
nisms, should the accelerated expansion in fact be due to 
the breakdown of Einstein gravity on cosmological scales. 
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